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The effects of locally broken spin symmetry are investigated in one dimensional nonequilibrium 
kinetic Ising systems via computer simulations and cluster mean field calculations. Besides a line 
of directed percolation transitions, a line of transitions belonging to N-component, two-offspring 
branching annihilating random-walk class (N-BARW2) is revealed in the phase diagram at zero 
branching rate. In this way a spin model for N-BARW2 transitions is proposed for the first time. 



I. INTRODUCTION 

The Ising model is a well known static, equilibrium 
model. Its dynamical generalizations, the kinetic Ising 
models, were originally intended to study relaxational 
processes near equilibrium states [fj]|2|. Glauber in- 
troduced the single spin-flip kinetic Ising model, while 
Kawasaki constructed a spin-exchange version for study- 
ing the case of conserved magnetization. Nonequilibrium 
kinetic Ising models, in which the steady state is pro- 
duced by kinetic processes in connection with heat baths 
at different temperatures, have been widely investigated 
and results have shown that various phase transitions are 
possible under nonequilibrium conditions, even in one di- 
mension (Id) (for a review see the article by Racz in Ref. 
H ) . Most of these studies, however, have been concerned 
with the effects the nonequilibrium nature of the dynam- 
ics might exert on phase transitions driven by tempera- 
ture. 

A different line of investigating nonequilibrium phase 
transitions has been via branching annihilating random 
walk (BARW) processes. The parity conservation of par- 
ticles is decisive in determining the universality class of 
the phase transition. A coherent picture of this scenario 
is provided from a renormalization point of view in [Q. 
The first example of a BARW model with an even number 
of offsprings exhibiting the so called PC (parity conserv- 
ing) transition was reported by Grassberger et al. 

A class of general nonequilibrium kinetic Ising models 
(NEKIM) with combined spin flip dynamics at T = 
and Kawasaki spin exchange dynamics at T = oo has 
been proposed by one of the authors Q in which, for 
a range of parameters of the model, PC-type transition 
takes place. This model has turned out to be very rich 
in several respects, for a review see 

Absorbing transitions have been, however mostly stud- 
ied in particle-type models. The N-BARW2 model is a 



classical stochastic system of N types of particles with 
branching annihilating random walk and two offsprings. 
For N — 1 the model exhibits, at finite branching rate 
p, PC type transition [8-14]. For N > 1 N types of par- 
ticles Ai perform diffusion, pairwise annihilation of the 
same species and branching Ai — > A{ + 2Aj with rate p 
for i = j and with rate p' /{N — 1) for i =/= j. In case of 
p = this model is always active except for the annihi- 
lation fixed point at zero branching rate. According to 
field theory the coarse-grained, bosonic version of the 
model forms a different universality class, the so called 
N-BARW2 with exponents in one dimension as follows: 
Ux = 1, z = 2, a = 1/2, P = 1 Here the exponents are 
defined as follows.: 



(i) 

(2) 



r is the characteristic time, £ is the correlation length 
and p(t), poo are the particle densities at time t and in 
the steady state, respectively. 

Hard core interactions have proven to be relevant in 
case of the N-BARW2 model by drastically changing the 
universality class |l7|]I§| ]. The arrangement of the off- 
springs relative to the parent turns out to be a relevant 
factor and causes two robust classes that are insensitive 
to parity conservation fl9|j or the binary nature of the 
production process p(J . 

In this paper we present an asymmetric spin-model 
(NEKIMA) with asymmetry both in the annihilation and 
spin-flip rate as a generalization of NEKIM. On the level 
of kinks, however, this model corresponds to a process of 
A and B particles with A — > ABA and B — > BAB -type 
branching and AB — > annihilation and BA — > BA 
exclusion. Nevertheless, as will be presented below us- 
ing computer simulations, the critical properties near the 
zero branching limit are the same as for the N-BARW2 
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model cited above with no sign of exclusion effects since 
alternating sequences of As and B-s occur like in [El], 
hence hard core interactions can not play an important 
role. Moreover, at finite branching rate of the kinks a line 
of directed percolation (DP) -type transition (l5) occurs 
which is well described by N = 6 level cluster mean-field 
calculations. 



II. THE MODEL 

The general form of the Glauber spin-flip transition 
rate in one-dimension for spin Sj sitting at site i is 
( Si = ±1): 



w(si, Si-i, s i+ i) = -(1 + 6si-is i+ i) 



(1 - -Sj(s;-i + s i+ i) 
(3) 



at zero temperature. Usually the Glauber model is un- 
derstood as the special case 6 = 0, T = 1. 

The Kawasaki spin-exchange transition rate of neigh- 
bouring spins Q at T — oo reduces to an unconditional 
nearest neighbour exchange: 



W ex (si, S i+ i) 



Pex 

2 



(1 - SjSj+l) 



(4) 



where p ex is the probability of spin exchange. The quan- 
tities in Eqs. ([?]) and (^) conserve spin symmetry, of 
course. Concerning spin exchanges, which act only at 
domain boundaries, the process of main importance here 
is that a kink can produce two offsprings at the next time 
step with probability 



Pk^3k oc p ex . 



(5) 



By changing p ex for negative values of 5 this model dis- 
plays phase transitions in the parity conserving univer- 
sality class [||. 

In the following we will be interested in investigating 
an extended version of the above model. Instead of (|^) 
we will prescribe the rates for the case T = 1, S = as fol- 
lows. The + + H > + + + and > processes 

remain as in (^) i.e. at zero temperature no kink-pair 
creations occur inside of domains. Further rates will be 
chosen in such a way that they break the symmetry of + 
and — spins locally. Such dynamically self-induced field 
was first investigated in a different context by Majumdar, 
Dean and Grassberger (MDG) Jig] , namely in studying 
the T = coarsening dynamics of an Ising chain in a 
local field which favors — spins as compared to + ones 
dynamically. 

In addition to the choice in Jl6[ | concerning the asym- 
metry in the annihilation rate 



w(+: ) = 1 

w(-;++) = 



(6) 
(7) 



further spin symmetry breaking will be introduced here, 
namely, in the spin-flip part of the Glauber transition 
rate the strength of which will be measured by a further 
parameter p + . While the transition rate 



«;(-,+-) = «;(-;-+) = 1/2 



(8) 



is unchanged, the two rates flipping + spins will be re- 
duced as 



w(+; + -) = w(+; -+) =p+< 1/2 



(9) 



in order to balance the effect of the other dynamically 
induced field arising from (6) and (7) by locally favoring 
+ spins. The spin-exchange part of the model remains 
as in the the spin-symmetric case, (Q). 

In the terminology of domain walls or particles the fol- 
lowing reaction-diffusion picture arises. There are two 
kinds of domain walls: — h = A and +— = B which can 
only occur alternately because of the spin background. 
Upon meeting AB — > while in the opposite sequence, 
BA, the two domain- walls are repulsive due to (7). 

The absorbing states in the extreme situation p + = 
when spin flipping maximally favors + spins, are states 

with single frozen — spins like H h + H hH h ++. 

By increasing p + , a slow random walk of these lonely — 
spins starts and by annihilating random walk only one 
of them survives and performs RW. The all + and all — 
states are, of course, also absorbing. 

The inclusion of nearest neighbor spin exchange, 
changes the picture drastically. Spin exchange leads to 
A — > ABA and B — > BAB type kink production, which, 
together with AB — * annihilation (BA — > is forbidden 
due to the annihilation-asymmetry, eqs. (6) and (7)) and 
diffusion of A and B leads to a kind of two-component, 
coupled branching and annihilating random walk. The 
phase diagram and the nature of the transitions will be 
reported and discussed in the following. 



III. PHASE DIAGRAM OF THE NEKIMA 
MODEL 

A. The line of DP transitions 

1. Simulation results 

Given three independent parameters 5, p+ and p ex (or 
rather p ex jT ) with the restriction eq. (B|) it is hard to 
explore the whole phase diagram. In the original (spin- 
symmetric) case of NEKIM [|| we have investigated the 
phase boundary in the parameter space (<5, p ex ). For neg- 
ative values of 5 a line of PC transitions was found. We 
are not going to investigate the S < case in the follow- 
ing and only make the remark that spin-asymmetry as 
introduced above (with a trivial generalization for 5 < 0), 
changes the parity conserving character of the transitions 
to directed percolation - type (of two species), as could 
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be expected. The case 5 > (including the Glauber case 
r = 1, S = Q) was found to be Ising-like, for all values of 

Pex • 

Introducing spin-asymmetry, however, makes the 
Glauber case also richer in phases. The phase diagram 
in the plane of parameters p ex and p + < 1/2 as obtained 
by computer simulations is shown on Fig.|l|. 




0.1 0.2 0.3 0.4 0.5 0.6 



P + 

FIG. 1. Phase diagram of the NEKIMA model for (5 = 0, 
r = 1. The absorbing (fully — ) phase lies above the bound- 
ary. The p+ — 0.5 point is referred to as MDG point in the 
text 

In fixing the phase boundary the quantity measured 
was the density of kinks, p(t), as a function of time start- 
ing from a random initial distribution of up and down 
spins. The chain size L varied between L = 2000 — 5000 
and up to t = 5 x 10 5 Monte Carlo steps (MCS) were 
reached. The way of updating was as in ref. §. A line 
of DP transitions has been found by the power-law be- 
haviour of pit) ~ t~ a with a — .160 ± .005, the value 
characteristic of DP transition. This was, of course, the 
expected kind of order-disorder transition on the basis of 
spin-asymmetry. 

The point p+ = 1/2, 

Pex — is a particular one, at 
this point our model goes over to the one investigated 
by Majumdar et al. pqj. For these parameter values a 
very high precision computer measurement has given the 
result p(t) = t^ 1 / 2 /ln(t/r) supporting an analytic inde- 
pendent interval approximation by these authors. In an 
equivalent model, however, the deviation from the i -1 / 2 
law has shown up as an initial density-dependent power 
function with a power of p(t) slightly deviating from 1/2 
|p2| . This problem, however, is not the subject of the 
present investigation. 

Concerning the phase diagram, Fig. (|l|), the param- 
eter values in the vicinity of p+ = 1/2 for p ex ^ are 
hard from the computational point of view as long tran- 
sients show up in the time evolution. It is apparent, 
however, that the phase line ends up here tangentially. 
At p+ = 1/2 the effect of the exchange term is such that 



for all p ex > the absorbing phase is entered: due to the 
choice in eqs.(6) and (7), the all Sj = — 1 phase (one of 
the absorbing phases) is reached exponentially fast. 

As to the other limiting situation p + — 0, for p ex = 
the initial spin distribution freezes in. As a matter of 
fact, the line of phase transitions reaches the p+ = axis 
only by letting p ex /^ — > oo by T — > ( T is only fixing the 
rate of flips, see eq.(3), while here we fixed it to unity). 
This circumstance, however, is of no importance for the 
results. 



2. Cluster mean-field calculations for the phase diagram 

Cluster mean-field approximation introduced for 
nonequilibrium models by |2^,Q was applied for the 
present model. The N — 1 mean-field equation for spin- 
up density is 

£gL = -2 P+ (l- Pl )p? (10) 

independently form p ex , which gives a p\ oc t~ x leading 
order decay to the pi(oo) = solution for all p + >= 0, 
while it is constant (keeps the initial value) for p + = 0. 
Therefore this predicts a discontinuous transition along 
the p + = axis. The corresponding steady state expo- 
nent is j3 = 0. Similarly the kink density, pi(l — pi), 
decays with a leading order singularity pkink oc t^ 1 and 
exhibits a jump at the pi n = axis. 

The N = 2 pair approximation results in the following 
steady state solution for kinks 



Pkink(oo) = 


4 Pex 


(l-2p+)p + 




1 + 2 (A Pex - 1) 


p+ + 8p ex (2p ex - 1 
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FIG. 2. Phase diagram determined by N = 2. ..6 cluster 
mean-field approximations (filled symbols from right to left), 
the iV — > oo extrapolated values (plus signs) and simulation 
results (stars). 
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This has absorbing states ( punk = 0) along the p+ = 0, 
p ex =0 andp + = 1/2 lines and active in the < p+ < 1/2, 
p ex > region. The transitions, however are continuous 
with leading order singularity (3 = 1 everywhere. As we 
can see the simple mean-field and higher order cluster 
mean-field approximations give different singular behav- 
ior similarly to cases treated earlier [ ^5|j2^j2^ , |29|] . 

For higher order approximations N > 2 we could 
solve the equations for the steady state numerically only. 
We could determine stable solutions up to the N = 6 
level from the coupled non-linear equations of 36 vari- 
ables. By locating the phase transition lines we found 
that the p+ = and the p ex = transitions do not 
change but the p+ — 1/2 (DP) transition line shifts 
monotonically towards the p+ = axis as we increase 
the level of approximations (see Fig. (|J). These so- 
lutions converge towards the phase transition line de- 
termined by simulations. We found that fairly good 
quadratic fitting can be applied for the N — 3, 4, 5, 6 level 
p* + {pex, N) critical point solutions, so we extrapolated to 
N — > oo at p ex — 0.1, 0.3, 0.5, 0.7, 0.9. The corresponding 
p* + {Pexi oo ) curve agrees well with the simulation data 



(see Table IV 



B. The line of N-BARW2 transitions 

As we have seen in the previous subsection, in the plane 
of [p +1 p ex ) the phase below the phase transition line is 
the active one and extends down to the p ex — axis. The 
critical behavior at and in the neighborhood of this axis 
has turned out to be of N-BARW2 type. In this respect 
the absorbing states are fully ordered or consist of single 
— es performing random walk in the sea of + spins. 
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FIG. 3. Effective critical exponent of the order parameter 
at the p ex = transition. 

The order parameter also in this case is the density 
of kinks p, the steady state value of it disappears when 



approaching the p ex — axis as p^ oc p@ x . Simulations 
from random initial state in a system with size L = 10 s 
were run up to 10 6 MCS. In the supercritical region the 
steady states have been determined for different p ex val- 
ues. Following level-off the densities were averaged over 
10 4 MCS and 1000 samples. By looking at the effective 
exponent defined as 



Peff(Pex(i)) 



In Poo(Pex(i)) ~ InPocjPexji - 1)) 
In p ex (i) - lnp ex (i - 1) 



(12) 



one can read-off: (3 e ff —> (3- The result of computer 
simulations at p+ = 0.1 is shown on Fig.(||). A linear 
extrapolation for p ex < 0.1 gives (3 = 1.0 ± .01. The 
overshooting of /? e // near the critical point is typical in 
case of logarithmic corrections to scaling. By plotting 
Poo I Pex as the function of la(p ex ) fairly good linear be- 
havior could be observed in the 0.02 < p ex < 0.4 region. 

The pit") simulation results at p ex = and p + = 0.1,0.4 
were analyzed by the local slopes 



*//(*) 



-\n[p(t)/p(t/m)} 
ln(m) 



(13) 



(where m = 8 is used). The asymptotic time evolution of 
the density of kinks p(t) ~ t~ a has proven to be, within 
error, that of annihilating random walk: a — 1/2 as 
shown on fig.(|4|). A ~ t~ 9 correction to scaling gave 
best fit in both cases. We also tried to fit a logarithmic 
correction form ((a+bln(t))/t) ' 5 for pit) but b was found 
to go zero for t >~ 6 x 10 5 MCS in both cases. 
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Effective critical exponent of the kink density at 
0, p+ = 0.1 (lower curve) and p+ — 0.4 (upper 
curve) as the function 1/t 0,4 . This choice results in linear 

' 9 correction 



FIG. 4. 

the pex = 



plot of the local slopes, corresponding to 
to scaling. 



t~ 



As to the remaining critical exponents when approach- 
ing p ex = we checked the expected N-BARW2 be- 
haviour by measuring the kink density in the active state 
(steady state) for several small values oip ex (between .1 
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and .001) for lattice size L between 50 and 5000. Th 
initial state was prepared in such a way that a cluste 
of '— '-es was chosen of width and location randomly dis 
tributed between L/4 and 3L/4. Finite size scaling the 
ory |3(| predicts the form 
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p(p ex ,L)=L-^F( Pex L 1 /^) 



(14 



Using the value of j3 as obtained above, we determin 
v± by data-collapsing. With f3/v± = 1 on Fig.(||), w 
get v± — 1.00 ± .06. Thus our result shows the critic;; 
exponent values by Cardy and Tauber 8. 
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IV. SUMMARY 

In this paper a new model has been presented by gen 
eralizing the nonequilibrium kinetic Ising model for th 
case when two kinds of spin anisotropics are present. In 
addition to the local kinetic bias introduced by Majum- 
dar et al [fl6| first, by prescribing with probability equal 
to zero the annihilation of — spins in the neighborhood 

'H h', we introduce spin-anisotropy in the spin-flip rate. 

Namely, the '+' spins are less likely to flip at domain 
boundaries (with probability p + ) than '— '-es (probabil- 
ity 1/2). Branching of kinks (domain boundaries) is the 
main effect of the spin-exchange part of the model. Wc 
have shown by cluster mean field calculations and com- 
puter simulations that for a given p + < 1/2 asymmetry 
the presence of spin exchange gives rise to two different 
types of phase transitions. While at p* x — an active 
phase emerges with an N-BARW2 type of transition, for 
p* ex > this transforms back into an absorbing state with 
a DP class transition. The spin anisotropies result in a 
new type of two-component, coupled branching and an- 
nihilating random walk of kinks with parity conservation. 

At p + = 1/2 the MDG point is reached, see Fig. 1. It 
is the endpoint of the DP transition line similarly to the 
compact directed percolation endpoint of the DP transi- 
tion line in the Domany-Kinzel cellular automaton model 
fjjj . The absorbing phase is the same in the two models, 
the active phase, however, is different. 

The question arises whether the 1 / ln(i) factor in the 
asymptotic behaviour found in |l6| ] is to be expected to 
hold even in the present model for p ex = 0. If this were 
the case, the N-BARW2 behavior would also be affected. 
First, from the side of simulations, we have not found any 
sign of such behavior. Moreover, the physical picture be- 
hind the expected asymptotic behaviour of spins is also 
different in the two cases. While MDG argue that at late 

time the process — I ► — with probability unity leads 

to '+' domains sandwiched between much larger '— ' ones 
the introduction of the local asymmetric spin-flip mag- 
netic field with bias for '+' spins will act against and feed 
up the '+' phase to compensate for their biased annihi- 
lation via the MDG process. As a consequence there is 
no reason to expect a late time logarithmic relaxation of 
the kink density on the p ex = line for p + < 1/2. 
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FIG. 5. Data collapse of pL 13 ^ 1 - against p e xL 1/v± - with 
f3/v± = 1 for various values of the chain- lengths (circle 
L = 50, diamond L — 100, plus sign L = 600, triangle 
L = 800, triangle left L = 1000, star L = 2000, triangle right 
L — 3000, triangle down L = 4000) on a double logarithmic 
scale. 

Finally it is worth noticing that while the PC transition 
is known to be sensitive to the symmetry |l4|]27j and 
DP transition appears by destroying it, the N-BARW2 
transition seems to be insensitive to this symmetry break- 
ing. 



TABLE I. Summary of p* + critical point results of N = 3—6 
GMF approximations and simulations 



pex 


N = 3 


N = 4 


N = 5 


N = 6 


N -> oo 


MC 


0.1 


0.419 


0.4003 


0.3903 


0.384 


0.357 


0.405 


0.3 


0.3357 


0.2963 


0.2772 


0.2661 


0.233 


0.272 


0.5 


0.3000 


0.2479 


0.2242 


0.2111 


0.178 


0.187 


0.7 


0.2824 


0.2216 


0.1945 


0.1798 


0.146 


0.137 


0.9 


0.2726 


0.2057 


0.1759 


0.1597 


0.122 


0.107 



TABLE II. Summary of critical exponent estimates at the 
Pex = line. The last row shows the data of the Id N-BARW2 
class. 



p+ 


P 


a 




0.1 


1.00(1) 


0.505(5) 


1.00(6) 


0.4 




0.503(5) 




N-BARW2 


1 


1/2 


1 
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